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On the coupling constant bounds implying area decay of Wilson loop for Z2 lattice 
gauge theories with Wilson action Boltzmann factor given by bond variables 

described by transverse Ising model. 
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We obtain correlation functions for lattice gauge theories with action Boltzmann factor given by 
bond variables described by transverse Ising model and apply them to obtain area decay of the 
Wilson loop observable in a range of the coupling parameter larger than that obtained from mean 
field theory considerations. 
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In this paper we study the lattice gauge theories using 
correlation identities and rigorous inequalities to obtain 
coupling parameters which imply area decay in lattices 
decorated by bond variables of spins a = 1/2 described 
by the transverse Ising action. Therefore we study the 
effect of a disordered field on the area decay behavior of 
the Wilson loop observable. Previously |l|, it has been 
presented the study of the quantum transverse Ising anti- 
ferromagnet on the kagome lattice, with weak transverse 
field dynamics and other local perturbations. The possi- 
bility of a disordered zero-temperature phase is investi- 
gated by means of an appropriate mapping to a compact 
U(l) gauge theory on the honeycomb lattice that is cou- 
pled to a charge- 1 matter field. The interest on such 
systems lies on the possibility of exotic phases and spin 
liquid states appearing in geometrically frustrated quan- 
tum magnets, as, for example, in Spin ice 
its mapping on a Ising system on a kagome lattice 
and excitations on the spin ice ground state [11 1 
constitute recent examples. In all these cases, in order 
to identify phase transitions on the system, the absence 
of a non-zero magnetization to be employed as an order 
parameter can be circumvented by evaluating the decay 
law of the Wilson loop [12j . 

Here, we consider the well known Wilson loop observ- 
able of a pure gauge Z d lattice theory with Wilson action 
defined as 



The action is, 



< W(C) >= lim < W(C) > A 



(1) 



where < W(C) >\, is the finite lattice Gibbs ensemble 
average, with action Boltzmann factor [l3T - [l5l ]. given by, 



(2) 



PCA 



where P denotes the unit squares (plaquettes) of A, and 
Xp is given by, xp = cr i cr l cr 3 cr l; °f an d °"f are the Pauli 
matrix for the spin 1/2. 
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Xp 



(3) 



PCA 



Let W(C) be the product of where b stands for 
the bond, along the perimeter of the planar rectangle 
C of area A. The gauge coupling constant is J , < 
j3J < oo. £1 is the local transverse field in the bond 
that exchanges its value from ±1 to =pl .Therefore the 
sign of the plaquette changes according to the ilaf term 
of the action Eq.((3|). Area decay of < W(C) > is a 
criterion for confinement. In section 3, we deduce the 
identity for < a z D >, where a z D = a^.-.a^ denotes a 
product of distinct bond variables. Similar procedure 
has been applied previously for the Z2 lattice gauge Ising 
case We take free boundary conditions and note 

that Griffith's first and second inequalities for the spin 
model apply [I8I - I20I ] and therefore imply the existence 
of the thermodynamic limit (2lj . For d = 3 and d = 
4 we obtain lower bounds (3l on the area decay of < 
W(C) > ( i.e. for aU /3 < fi L , < W(C) > has area 
decay, using correlation equalities for <j z D and Griffith's 
inequalities. For completeness, in section 2, we give the 
mean field lower bound /3m (Pm < Pl ) using a decoupling 
and Griffith's inequality argument [16l |22j . 

For definiteness assume C lies in the x\ , X2 plane. Con- 
sider a bond b fixed in the lower left-hand corner of C. 
Replace (3J by ((3J)\, A £ [0,1], in the action for the 
2(d—l) plaquettes (call them Pi.-P^td-i)) that have one 
bond in common with b. In the following we will discard 
the mechanism of disorder represented by Q. Denote the 
corresponding expectation by < W(C) >\ . Integrating 
(-£■) < W{C) > x gives, noting that < W(C) > = 0, 



f 1 d 

< W(C) >= / d\— < W{C) >> 
Jo d\ 

f l 2(d-l) 

= dX V (< W(C) XPi >x 
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< W(C) > A >< X.Pi >> 



2 



2(d-l) 



< p [ rfA V < W(C) XPl >) 
Jo i=i 



2(d-l) 



< p £ < W(C) XPi >, 



(4) 



using Griffith's first (second) inequality in the first (sec- 
ond) equality. Each term on the right corresponds to a 
modified contour determined by the bonds of the vari- 
ables of W(C)xPi which enlarges or diminishes C by one 
plaquette. We repeat the argument proceeding along suc- 
cessive rows of plaquettes enclosed by C. After A appli- 
cations we arrive at 

< W{C) >< /3 A (sum of [2(d - 1)] A terms) (5) 



we obtain, writing H = H b + H' , where H b includes all 
terms related to bond b and H 1 the rest of the lattice, 
the equation, 



z{b)Tr (b) <r z D e- pHb 



D 



Tr 



(b)< 

Tr (b) a z D e'^ 



-0H b 



Tr 



(b)P D t 



(11) 



where Tr/. 



-(h) represents the partial trace with respect 
to bond b and A = 1 - e -PH be -pH' e p(H b +H') _ Equation 
(fTTj) is an exact relation. However, it is difficult to be 
used. Therefore, we will make an approximation based 
on the following decoupling, 



Each term is non-negative and bounded above by 1 giving 
W{C)\ < [/32(d - l)] A . Therefore, for each /3 , 



we have, 



/3ei/(2(d-i)), 



W(C)^ < e [-log(2^(d-l)A)] > 



(6) 



(7) 



In the preceding arguments, which led to the mean field 
lower bound Pm, we have modified the action by replac- 
ing pj by (/3J)A. The effect of O is to lower the criti- 
cal temperature T c , or to increase j3 c , up to T c = 0, or 

Pc = OO. 

We now present the derivation of the longitudinal spin 
correlation function which is derived in a manner anal- 
ogous to the ones for the spin transverse Ising systems 
[23| (see also [24j) and they will be used to extend the 
P region of area decay given by relation ([5]), for d = 3 
and 4. Let of, = u z x ...a z D denote a product of distinct 
bond variables and for a fixed bond b occurring in a z D 
give a numerical ordering 1,2,... to the 2(d- 1) plaquettes 

that have one bond in common with b. Let a$ be the 
product <td with the bond b deleted. Let us consider the 
bond b, with a b , and the plaquette Xb which contains the 
bond b. Summing over a b , we have, 



< <r D > = 



(t-q 1 — — — tanh pEi 



E, 



where, 



with Ej 



Ei = JE. 



E 2 



(8) 



(9) 



and Ei 



Also, tjj,crt 



and of are the neighbours of of on the 



plaquette Xb- 

From the expectation value 

< of, >= Tr(o* D t 



(10) 



Mb) 



z(b) 



Tr {h) o z D e-^ 



Tr 



Tr [b) <7 z D e-^ 



Tr 



(12) 



Inserting (fT2|) into (JTTJ) and using the fact that A ^ 1. 
we obtain, 



< On >< 



z(b) Tr(b)V z L 



-PH b 



Tr 



(&)< 



-PHb 



(13) 



By expanding A we see that the approximation is cor- 
rect to the order of p 2 and it is consistent with the ap- 
plication of the correlation inequalities that will be used 
later. Remarck : A more general result can be repre- 
sented by equations identical to (ITT1) and (flU)) but writ- 
ten for < F(a)a b >, where F(o) is any product func- 
tion of the bond components, except o|. Relations (fTT|) 

and (flU)) are special cases where F(a) = o^, . The re- 
lation for the longitudinal bond magnetization < o^ > 
is the special case where F(a) = 1. Within this de- 
scription, the decoupling (fl~2"|) can then be viewed as a 
0-th order approximation of the exact relation (1111) for 
F(S) = 1 , i.e. , Eq. (fT3"|) can be assumed to be ob- 
tained from the approximation < o^ — o^ >= where 
^ b =Tr (b) {a z b e-^}/Tr (b) {e-^}. 

Introducing V 
get, 



dx 



through e aV ./(x) = f(x + a), we 



< °£> >< 



fix) U4<J4) 



and, 
fix 



vm 2 

Considering, 



: tanh(^) V / (2fi) 2 + x 2 = -f(-x) 

(15) 



3 



' i,k,l 



then (W(C)) < 



for d= 3. 



(21) 



j,k,l 
we obtain, 



Y[ [cosh(J 6j -ikiV) + o]alal eunh(J ww V) )jl6) 



< ah >< Yl [cosh(Jy H V) 

+<r^af smh(J 6iW V)] )./(*) U =0 (17) 



Applying the previous result, equation (|17l) . for d= 
2,3,4, after some algebra, we obtain, (a), (b) and (c): 

(a) d = 2 

< <rf, >< a 2 < CT c^ >' ( 18 ) 



a 2 = -/(2J),a 2 >0 



(6) d = 3 

< erf, >< a 3 J2 < a hXPi > 

i 

+b 3 Y, < vhXPiXPjXP* >, (19) 

ijk 

«3 = ^[/(4J) + 2/(2J),a 3 >0 
6 3 = ^[/(4J)-2/(2J),6 3 <0 



(c) d = 4 

< erf, >< a 4 ^ < cr^XPi > 

i 

+b 4 ^ < vhXPiXPjXP* > 

+c 4 ^2 < ahXPiXPjXPkXPiXPm >> (2°) 

ijklm 

«4 - ^[/(6J) + 4/ (4 J) + 5/(2 J), a 4 > 

^4-^[/(6J)-3/(2J),6 4 <0 

c 4 - ^[/(6J) - 4/(4 J) + 5/(2 J), c 4 > 

We will now derive the coupling bounds for the area 
decay of the correlation functions . Applying the proce- 
dure described below to equations (I2U1) and (|2"Tj) . we will 
obtain the upper bounds for d=3 and d=4. If ft J is such 
that: 

(a)4a 3 < 1, 



(6)6(a 4 + c 4 ) < 1, 

then /W(C)) < e (-'«[ 6 (^+c4]vi) )for d= 4 ( 2 2) 

The proofs of results (a) and (b) follows the same pro- 
cedure used to obtain the mean field results (section2). 
Let us prove result (a). For b e C, 

< W(C) ><a 3 Y,< W ( G )XPi > 

i 

+b 3 J2 <W{C) XPt XP 3 XP k > 

i<j<k 

<a 3 Y / <W(C) XPt > (23) 

i 

The last inequality is obtained because b 3 is negative 
and < W (C)xPiXPjXPk > i s positive by Griffith's first 
inequality. At each application of the equality we pass 
to an inequality by dropping the 63 terms. After A steps 
we arrive at 



< W{C) >< af{ sum of 4 terms) 



(24) 



where each term is less than one. To prove result (b) for 
d=4 we use equation (|20[) and proceed as previously as 
in the proof of result (a) for d=3 . We can drop the 6 4 
terms in favour of an inequality. At each stage we have 
six terms from the a 4 term and six terms from the c 4 
term. 

The values of £l c , the local disorder critical field, for 
which the conditions (|22p and (f2"2"j) are not satisfied, in 
other words, for which the perimeter decay regime of the 
Wilson loop dominates, are: 



n c (d : 
fl c (d : 



3) = 1.3755J 

4) = 2.4466J 



(25) 
(26) 



For values greater than those there are no confinement. 

In this paper, we introduce the Z d lattice gauge model 
for which the action has a term that represents a disor- 
der mechanism. This term is described by a local field 
represented by the a x spin 1/2 operator. The effect of 
this local field is to change the sign of the plaquette, 
i.e, from ±1 to qpl. We have use correlation identities 
and rigorous inequalities, to obtain, for d = 3 and d = 4, 
lower bounds for the area decay of < W(C) >, i.e , for all 
ft < fth, < W(C) > has area decay. The correlation iden- 
tities are a gauge version of Callen's identities employed 
by the authors [23j to obtain lower than mean field upper 
bounds for the critical couplings for the transverse Ising 
model. The procedure was based on an approximation 
for an exact identity and on rigorous inequalities for the 
spin correlation functions. We obtain the critical local 
disorder field above which the perimeter decay regime of 
the Wilson loop dominates. 
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